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Factorization of 

Multivariate Positive Laurent Polynomials 

Jeffrey S. Geronimo^^ and Ming-Jun Lai^^ 


Abstract. Recently Dritschel proves that any positive multivariate Laurent 
polynomial can be factorized into a sum of square magnitudes of polynomials. 
We first give another proof of the Dritschel theorem. Our proof Is based on 
the univariate matrix Fejer-Riesz theorem. Then we discuss a computational 
method to find approximates of polynomial matrix factorization. Some numer¬ 
ical examples will be shown. Finally we discuss how to compute nonnegative 
Laurent polynomial factorizations in the multivariate setting. 


§1. Introduction 

We are interested in compnting factorizations of nonnegative Lanrent polynomials 
into snm of sqnares of polynomials. That is, let 

n 

p{z) = 

k=—n 

be a Lanrent polynomial, where 2 ; = Suppose that P{z) > 0 for \z\ = 1. One 

n 

would ask if there exists a polynomial Q{z) = qkZ^ such that 

fc =0 

P{z)=Q{zrQ{z), ( 1 ) 

where Q{z)* denotes the complex conjugate of Q{z). This is the well-known Fejer- 
Riesz factorization problem and it was resolved by Fejer [F’15] and by Riesz [R’15]. 
A natural question is whether the results of Fejer and Riesz can be extended to the 
multivariate setting. More generally, given a nonnegative multivariate trigonomet¬ 
ric polynomial P{z) := P{zi, Z 2 , - ■ ■, Zd) of coordinate degrees < n, does there exist 
a hnite number of polynomials Qk (z) such that 

= ( 2 ) 

k 
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i.e., can P{z) be written as a snm of sqnare magnitndes (sosm) of polynomials. 
There is a vast amonnt of literatnre related to the stndy of this problem and the 
resnlts relevant to this paper may be snmmarized as follows: 

1° When P{z) is nonnegative on the mnlti-torns l^il = |^ 2 | = ■ ■ ■ = \zd\ = 1 and 
the coordinate degrees of Qk are less than or eqnal to n, then the answer to 
the qnestion is negative. (See [Calderon and Pepinsky’52] and [Rndin’63].) 

2° When P{z) is strictly positive on the mnlti-torns and the coordinate degrees 
of Qk are not specihed, Dritschel has shown that the answer to the qnestion is 
positive([Dritschel’04]). However the nonnegative case remains nnresolved. 

3° In the bivariate setting, Geronimo and Woerdeman gave a necessary and snfh- 
cient condition in order for P{z) = \Q{z)\‘^ ^ where Q{z) is a stable polynomial, 
i.e., Q{z) 7 ^ 0 inside and on the bi-torns ([Geronimo and Woerdeman’04]). 

4° In the bivariate setting, there exist rational Lanrent polynomials Qk{z) snch 
that (2) holds. Fnrthermore, Qk can be so chosen that the determinants of Qk 
containing only one variable Lanrent polynomials (cf. [Basn’Ol]). 

5° In [McLean and Woerdeman’Ol], an algorithm was proposed to hnd polynomi¬ 
als Pk snch that P = J2k algorithm nses the so-called semi-dehnite 

programming. 

Althongh the mathematical problem appears to be theoretical, it has many 
applications in engineering, e.g., the design of antoregressive Liters, constrnction of 
orthonormal wavelets (cf. [Danbechies’92]), constrnction of tight wavelet framelets 
(cf. [Lai and Stoeckler’04]), spectral estimation in control theory (cf. [Sayed and 
Kailath’Ol]) and many other engineering applications mentioned in [McLean and 
Woerdeman’Ol]. Thns, how to compnte snch factorization polynomials Qi, Q 2 , ■ • •, 
is interesting and nsefnl for applications. In this paper, we discnss a symobl approx¬ 
imation method stndied in [Lai’94] for compnting snch factorizations. The method 
was originally intended for factorizing any nonnegative Lanrent polynomials in the 
nnivariate setting. We nse the ideas to give Dritschel’s theorem another proof. The 
proof provides a compntational method to factor P{z) into Qk{zys. The paper is 
organized as follows. In section 2, we hrst give a different proof of Dritschel’s The¬ 
orem. A key in the proof is to factorize nnivariate Lanrent polynomial matrices. In 
section 3, we discnss how to compnte the factorization of positive Lanrent polyno¬ 
mials matrices. The method nsed by [Lai’94] to compnte approximate factorizations 
is extended to the matrix case. Then in section 4, some n um erical examples are 
compnted following the procednre in §2 and §3. Finally in section 5. the nonegative 
case is considered. 


§2. Dritshel’s Theorem 

We begin with reviewing the concept of the symbols of bi-inhnite Toeplitz matrices 

n 

and its properties. For a given nnivariate Lanrent polynomial P{z) = E Pkz'", 

k=—n 

we may view P{z) as the symbol of a bi-inhnite Toeplitz matrix V := {pi-j)ij^z- 
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Indeed, for any absolutely summable sequence x = {xi)i^z, be., \xi\ < oo, let 

i^Z 

F(x) = discrete Fourier transform (or z-transform) of x. Let 

y = "Px, then it is easy to see that 

F(y) = P(z)F(x). 

If the matrix V has a factorization Q which is a banded upper triangular Toeplitz 
matrix such that 

'P = Q^Q, 

the discrete Fourier transform of y = QQ^^x is -F'(y) = Q(z)*Q(z)F(x), where Q) 
denotes the complex conjugate transpose of Q. Thus, hnding P{z) = Q{zYQ{z) 
is equivalent to hnding a banded upper triangular Toeplitz matrix Q such that 
V = Q^Q. 

It is easy to show that if P{z) > 0 for all \z\ = 1, then V is Hermitian and 
nonnegative dehnite. Clearly, V is Hermitian since P{z) is real. Furthermore for 
any absolutely summable sequence x, we need to show that x^Px > 0. Again 
writting y = Px, we know that 

1 /‘ 27 '-_ 

xV = ^ P{x)P{y)de 
where z and it follows that 

xtpx=— / \P{x)\‘^P{z)d9 >0 

27r Jo 

for any nonzero sequence x. In particular, for 

X = (•■•, 0, X-N, • • • , Xo, • • • , XJV, 0, • ■ 

the left-hand side in the above inequality gives x^P/vx, where P/v is a central section 
of P. The above argument shows that Pn is nonnegative dehnite. 

In the following we will assume that P{z) is strictly positive, in the sense that 
there exists a positive number e > 0 such that P{z) > e. When P{z) is a matrix, 
we mean that P{z) > el, where I is the identity matrix of the same size as that of 
P{z). When P{z) is strictly positive, we have 

xtpx = — / |P(x)|2p(^)d6' > e||xf. 

27r Jo 

It follows that if P{z) > e > 0, then P/v > e > 0. 

We now consider the factorization of multivariate Laurent polynomials. Let us 
begin with a bivariate Laurent polynomial P{zi, Z2) hrst. That is, let 

n n 

P{zi,Z2) = > 0 

j=—nk=—n 
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be a Laurent polynomial of coordinate degrees < n. We would like to find a finite 
number of polynomials Qk such that 

P{Z1,Z2) = \Qk{zi,Z2)\‘^. 

k 


Denote by zi = [1, zi, zf, • • •, Zi]'^ and write 


P{zi,Z2) = Z1^P(Z2)Z1 

n 

for a Hermitian matrix P{z 2 ) = PkZ^^ where each pk is an (n + 1) x (n + 1) 

k= — n 

Toeplitz matrix. With a slight modihcation of an observation of [McLean and 
Woerdeman’Ol, Theorem 2.1], we note that there are many ways to write P{z 2 )- 
If there is one P{z 2 ) which is nonnegative dehnite then we can use the matrix 
Fejer-Riesz factorization (cf. e.g., in [Helson’64], [Mclean-Woerdeman’Ol], see also 
section 3) to hnd < 5 (^ 2 ) such that 

P{Z2) = QH^2)Q{z2)- 


That is, we have 

P{zi,Z2) ^ iQ{z2)ziyQ{z2)z:L 

which is clearly a sum of squares of polynomials. 

The above discussion can be generalized to the multivariate setting and using 
an observation of [Dritschel’04] to the case that the size of P{z 2 ) is larger than 
(n + 1) X (n + 1). For simplicity, let us consider a trivariate Laurent polynomial 
P{zi, Z 2 , zs) in zi = ^2 = ^ zs — of coordinate degrees < n. We hrst 

write P{zi, Z 2 , Zs) in a matrix format: 


n 

P{zi, Z2, Zs) = YjPk{z 2 , Zs)z'l = Zi^P(z2, ^ 3 )zi, 

—n 


with 

= (3) 

and mi > n. There are many ways to write P{z 2 , zs). To capture this dehne the 
set of matrices 

P{z 2:Z3) ^ {{pij{z2,Z3)) 0 <ij <mi : Y Pi,jiz2:Z3) ^ Pk{z2,Z3),}. 

i-j = k 
\k\ < mi 
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Note that the matrices in T are banded since pk = 0, \k\ > n. We look for a matrix 
P{z 2 ,Z 3 ) in T that is positive dehnite for \z 2 \ = 1 = \z 3 \. The polynomial matrix 
P{z 2 ^ Z 3 ) can be written as 


n 

P{Z2,Z3) = ^ Pk{z3)z2, 

k= — n 

where each Pkizs) is an (mi + 1) x (mi + 1) Toeplitz matrix. Thns we can write 


^(^ 2 ,^ 3 ) = Z2^P{Z3)Z2: 


where 

with Im^ being the (mi + 1) x (mi + 1) identity matrix and m 2 > n. The polynomial 
P{z 3 ) is a matrix polynomial of size (mi + l)(m 2 + 1 ) x (mi + l)(m 2 + 1). If 
it is nonnegative dehnite we can factor it into a polynomial matrix Q{z 2 ), i.e., 
P{z 3 ) = Q{z 3 y Q{z 3 ) by the matrix Fejer-Riesz theorem (cf. [Helson’64] or [Mclean 
and Woerdeman’Ol]) then we have 


P{zi,Z2,Z3) = (Q(^3)z2Zi)^ (Q(^3)z2Zi) 

which is a snm of sqnare magnitndes of polynomials in zi, 2 : 2 , ^ 3 . 

Onr task then is to prodnce a positive dehnite polynomial matrix for any given 
positive mnltivariate Lament polynomial. We resnme onr discnssion on the two 
variable case again and rewrite P{zi,Z 2 ) as follows: 


ni 

-P(^1,^2)= ^ PA:(^2)^^ = Zmi^Tbii(^2)Zmi 

k=—ni 

where mi > ni, Zmi = [I7 ‘ ‘ 7 ^ and 


Pm I (^ 2 ) \Pjk{z2^]o<j,k<mi 


with polynomial entries Pj,k{z 2 ) given by 

Pjk{z2) = - , , ^ I . - ^Pk-j{z2),yj,k = 0,---,mi. 

mi + 1 — \j — k\ 

Note that Pjk{z 2 ) — 0 for \j — k\ > ni. Under this decomposition we can show 
that for some mi large enongh, the matrix Pi ( 22 ) will be positive dehnite when 
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P{zijZ 2 ) is positive definite. To see this we note P{zi,Z 2 ) is the symbol of the 
following bi-infinite Toeplitz matrix, 


P0{Z2) 

P-1{Z2) 

Pl{z2) 

P0{Z2) 

P2{z2) 

Pl{z2) 

Pn{z2) 

Pn-liz2) 


P-n{z 2 ) 0 


(4) 


The positivity of P{zi,Z 2 ) implies that any central section of the this matrix, i.e., 
any sqnare block with the diagonal consistent with the main diagonal 

diag(- ■ ■,po{z 2 ),Po{z 2 ),Po{w 2 ), ■ ■ ■) 


is positive as explained at the beginning of this section. Typically, we have 


Po{z 2 ) > 0, 


Poiz2) P-liz2) 
Pliz 2 ) P 0 {Z 2 ) 


P0(Z2) P-l(z2) P-2 (z2) 
Pl(z2) Po(z2) P-i(z2) 
P2(z 2) Pl(z2) P0(Z2) 




For convenience, we denote by P 2 and P 3 to be the 2x2 and 3x3 matrices above, 
respectively. In general, we nse Pk to denote the k x k central block matrix from 
the bi-infinite Toeplitz matrix (4) above. 

We now look at the matrix Pmi{z 2 ) given by. 


1 


mx + l 


P 0 {Z 2 ) 




m\ 


mi 


^^ 2 (^ 2 ) 


T^P 0 {Z 2 ) ^^P- 1 {Z 2 ) 


1 

mi 


Ti(^2) 


(^ 2 ) 

0 

: ^Poi^ 2 )- 

Note that each diagonal snms to Pi-j{z 2 ) so Pm^ is in P{z 2 ) where P{z 2 ) is defined 
as above with the obvions modifications. With x = [xq, xi, • • •, 0 :^, 1 ]^, we need to 
prove that x*Pm^{z 2 )x > 0. First we write 


X 


^PmA^ 2 )x = 


1 


mi -|- 1 


-X' 


-I- 


1 


mi -|- 1 


x' Rm^x 
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with a remainder matrix Rmx- The £2 norm of Rm^ can be estimated using the 
column norm to give 

ni(ni + l)Ci 

^ “ 2 (mi — ni) 


where Ci = sup^ |pi(z2)|- If P{zi^Z2) > e then x'^VmiX > e||a:||2, so that 

if ^ < e, then (^2)3:^ > 0 . Then an application of the matrix 

2(mi — ni) 

Riesz-Fejer Theorem yields 


ni 

Theorem 2.1. Let P{zi,Z 2 ) = Pfc(^ 2 )^i > e > 0 be strictly positive on 

k= — ni 

bi-torus \z\ = 1 — \z 2 \- Then P(zi,Z 2 ) can be factored into a sum of squares of 
polynomials in zi and Z 2 - The total number of terms in the sum is less than or 
equal to mi + 1 with mi being an integer such that 


ni(ni + l)^! 

2 (mi — ni) 


and the degrees of each of the polynomials is bounded by mi in zi and 77-2 in Z2 ■ 

We remark that when P(zi, Z2) has different coordinate degrees ni, 77-2, it may 
be worthwhile depending upon Ci to choose the smaller among 771 and 772 in order 
to have a fewer terms in the sum of square magnitudes of polynomials for P{zi, ^2)- 

Next we generalize the result in Theorem 2.1 to the multivariate setting which 
is known from [Dritschel’04]. 

Theorem 2.2(Dritschel). Let P{zi,---,Zd) be a multivariate Laurent polyno¬ 
mial which is strictly positive on the multivariate torus \zi\ — \z2 \ — ■ ■ ■ — \zd\ — 1, 
where d > 2 is an integer. Then P{z, w) can be expressed as a sum of square 
magnitudes of polynomials in zi, ■ ■ ■, Zd- 

Proof: We shall use the arguments in the proof of the previous Theorem. Write 
P{zij Z 2 t ■ ■ ■ T Zd) = P{zi,z) = ^ where z is the usual multi- 

variable notation beginning with Z2. We know that P{zi,z) is the symbol of the 
bi-in£nite Toeplitz matrix given by (4) with Z2 replaced by the multivariable 2 . It 
follows that any central section along the main diagonal is strictly positive dehnite 
as explained before. Write 


P{zi,z) = P^^{z)zi, 


( 5 ) 


where zi given by equation ( 3 ) and Pmi{z) = \Pj,k]o<:j k<m ^ matrix of size 
(7771 -|- 1 ) X (7771 -|- 1 ) with entries 


Pjk 


-—-p- -Pj-kiz), Vj,/C = 0 , 1 , •••, 7771 . 

777 i -I- 1 — |j — /c| 
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If p > e the argument in Theorem 2.1 shows that for mi large enough there is an 
ei > 0 such that {z)x > ei on the d-1 torus if ^ ^< e, where in 

2mi — ni 


n2 


this case Ci = sup^ |Pi(^)l- Write P^^{z2,z') = E Pk{z')z2, where 

k= — n2 

pk are (mi + 1) x (mi + 1) Toeplitz matrices and z' — {zs, ..., Zd). Now set 


Pjk 


W2 + 1 - |j - k\ 


Pj-k{z'): Vj,/c = 0, ■••,m 2 


with 1712 > xni and Pm 2 {^') = [Pj,fc]o<j k<m 2 ' above we have that 


1 1 

rp 1" fy rp rp 1" ' rp I rp 1 "' rp 

^ mo'^ , -t ' TUo'^ I , -t ^ 

m2 + 1 m2 + 1 


As above the norm of 7lrn2 can be bounded by 




^2(^2 + 1)^*2 
2(m2 - 172 ) 


where C 2 = sup ||pi(^0ll2- Thus for m 2 sufficiently large, P^j is a positive 

i,\z2\ = ---\zd\ = 7 

matrix polynomial. We continue the process until we arrive at the positive trigno- 
metric matix polynomial Pm^.i (^d) which can be factored by the matrix Fejer-Reisz 
Theorem. We have thus established the proof. □ 

Note that the number of factors will be (mi + l)(m 2 + 1) ■ ■ ■ (m^-i + 1) and 
the degrees of the polynomials at most mi for 2:1 ... m^-i for Zd-i and Ud for Zd 
. We note that we could have avoided the use of the matrix Fejer-Riesz lemma by 
eliminating all variables then using a square root of a positive matrix (see Mclean 
and Woerdeman’Ol). We will consider an alternative computationally attractive 
method for computing factorizations in the next section. 

§3. Computing Approximate Factorizations 

As shown in the previous section, an important step in the factorization of multivari¬ 
ate Laurent polynomials is to compute the factorization of univariate polynomial 
matrices. Recall a computational algorithm for factorizations of one variable trig- 
nometric polynomials was exploited in [Lai’94]. This method can be developed to 
factorize polynomial matrices in the univariate setting. Let us hrst introduce some 
necessary notation and dehnitions in order to explain the method in more detail. 

Let £2 stand for the space of all square summable sequences. Let ||x ||2 denote 
the standard norm on £ 2 . We note that for any operator A from 1 —> A can be 
expressed by a bi-inhnite matrix. 



Definition 3.1. A bi-inGnite matrix A = (aik)i,kez is said to be of exponential 
decay off its diagonal if 

IKkh < 

for some constant K and r G (0,1); where Z is the collection of all integers. A is 
banded with band width b if aik = 0 for all i,k E Z with \i — k\ > b. 

If A is a positive operator, then there exists the unique positive bi-infinite 
square root matrix Q of A such that = A. li A = B for another bi-infinite 
matrix B, then there exists a unitary matrix U such that B = UQ. 

Recall from the previous section that given any Laurent polynomial P{z), we 
can view P{z) to be the symbol of a bi-infinite Toeplitz matrix V. The computa¬ 
tional scheme introduced in [Lai’94] roughly speaking is to choose a central section 

Pn = {Pj-k)-N<j,k<N 

of matrix V and compute a Cholesky factorization i.e Pn = C^^p^Cn where Cn is 
an upper triangular matrix with positive diagonal entries, if Pn is positive definite 
or use the singular value decomposition (SVD) to find Qat if Pn is nonnegative 
definite and then find a Householder matrix Hn such that Cn = HnQn is upper 
triangular. Then the nonzero entries in the middle row oi Cn approximate that in 
the middle row (in fact any row) of C whose symbol C{z) is a factorization of P{z), 
i.e., P{z) =C{z)*C{z). 

For the extension of this method to matrix polynomials, let 

= {x = {rrjiez, Xi E R"", ||x||fc < cx)}, /c = 1, 2 

and B{11£) be the set of bounded linear operators on . Let H^r G B{P^) be the 
projection given by 

HatX = y, y = {yi} : yi = 0, |f| > N,yi^ Xi, |f| < N. 

If P G B{i‘^) is positive definite we will be interested in considering the {2N + 
l)m X (2iV -|- l)m submatrix of P centered at the index zero which will be called 
the Nth central section and which is also positive definite. We will also be interested 
in extensions of various finite matrices An to B{P^) given by 

"O 0 o' 

0 An 0 , 

0 0 0 

which with a slight abuse of notation will also be called A^;-. 

Consider the matrix polynomial P{z) = with matrix coefficients 

Pfc’s of size m X m, then V = {pi-j)ij^z G B{fT) defined by m x m matrix blocks 
Pk: —n < k < n IS a bi-infinite block Toeplitz matrix whose symbol is P{z). As 
shown earlier if P{z) is Hermitian nonnegative definite, so is V. Let C{z) be a 
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factorization of P{z) i.e., P{z) = C'(z)^ C(z), then V = where C is a bi-inhnite 
npper triangnlar banded block Toeplitz matrix associated with C{z). On the other 
hand, if P = for a npper triangnlar banded block Toeplitz matrix, then the 
symbol C{z) of C satishes P{z) = C{zyC{z). If P{z) is positive dehnite then it 
follows from the matrix Fejer-Riesz Lemma [Helson’64], [Mclean-Woerdeman’Ol] 
that it is possible to choose C so that it has positive diagonal entries. We shall 
prove the following. 

Theorem 3.1. Let P{z) = 'Yl'^nP^z^ be an m x m matrix polynomial that is 
positive dehnite for \z\ — 1. Let V = where C is an upper 

triangular banded block Toeplitz with positive diagonal entries, Pn be the Nth 
central section of V, and Cn its Cholesky factor (which we extend as described 
above). Then 

\\iCN-CN)5h<Kp^, 

for some p G (0,1), where 5 & is a vector with a hnite number of nonzero entries. 

Remark For the n um erical compntation below we will choose S with zero compo¬ 
nents except for So = fm, the m x m identity matrix. 

The proof of Theorem 3.1 is based npon the following Theorem 3.2 and Lemmas 
3.3 and 3.4. 

Theorem 3.2. Suppose that A G 5(^2) is a positive banded operator such that 
||A — /II 2 < 1. Let Q be the unique positive square root of A, Ajy be a central 
section, and Qn be the positive matrix such that — Aj^. Then 

\\iQ-QN)S\\2<KX^ (6) 


for some A G (0,1) and a positive constant K. In equation ( 6 ) S E I2 is any vector 
with a hxed number of nonzero entries. 

To prove the above Theorem, we begin with the following lemmas. 

Lemma 3 . 3 . Suppose that A is banded with bandwith b and || A — /||2 < r < 1. 

Then Q = {qik)i,kez then \qiy\ < Kr~^. If A is invertible, then the entires of 
Q~^ satisfy a similar bound. 

Proof: We only prove the exponential decay property of Q. The proof of that of 
Q~^ is similar. The nniqneness of Q and the convergence of the following series 


i=0 


(2i-3)!! 


{A-iy 


implies that 


e = VI = /TTcavt) = - ly. 

i=0 
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A is banded and so is A — If A — I has bandwidth 6, then {A — ly is also banded 
with bandwidth ib. Thns, 


where {A — I)jk denotes the (j, ky^ entry of A — / and similar for {A — lyj^- It 
follows that 

\qjk\ < 

for some constant K. This hnishes the proof. □ 

Let ns write 



C(i 

B 

Ct2 


'Pi 

a 

P2 

Q = 

St 

Qn 

(7t 

and A — 

at 

An 

ct 


Cts 

C 

a4 


Ps 

c 

Pr 


Note that = A implies A^ = Qj^ + + C^C or Qj^ — Qj^ = B^B + C^C 

where = A^. Thns, we have 

(Qn + Qn){Qn — Qn) = On ~ Q% + QnQn — QnQn = + (7^(7 + R, (7) 

where R is dehned in the following, 

Lemma 3.4. (cf. [Lai’94]) Let R = {rjk)-N<j,k<N ■= QnQn — QnQn- Then 
rjk = for k = -N/4 + 1, • ■ ■, N/4 - 1 and j = -iV, ■ ■ ■ , iV. 

Proof: (of Theorem 3.2.) From eqnation (6) we hnd that, {Qn — Qn) = {Qn + 
Qn)~^{B^B + C^C + R). By Lemma 3.3., we can prove that the entries of + 
C^C have the exponential decay property: {B^B + C^C)jk = —N < k < 

N. 

The positivity of A implies that Q is positive and so is Qat- It follows that 
||(5)(f^||2 is nniformly bonnded, fnrthermore since hatQN ia also nonegative we hnd, 

IK^at + Qn) < IIQ^/'^lb < i7i < cxD 

for a positive constant Ki independent of N, where we have nsed the fact that Qn 
is nonnegative. Therefore, we conclnde that 

IKQiV - Qn)Sn\\2 <||(QAr + QN)-y\\\{B^B + C^C + i?)hAr||2 

<Ki\\{B^B + C^C + R)5Nh 

where Sn is the hnite vector whose entries match those of 6. The proof is completed 
by extending Qn, Qn, replacing Sn by h, and noticing that by Lemma 3.3 IKQat — 
g)h||2 < i 7 iA^, A < 1 . □ 
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Proof: (of Theorem 3.1) Suppose that 


sup ||P (^)||2 < 1. (8) 

1^1=1 

Otherwise divide P by a sufficiently large constant so that (8) holds. Let Q be the 
unique positive square root of P, and Qn the positive square root of Pn- From 
Theorem 3.2 we know that || (Qat — Q)h ||2 < Kp^ with p < 1. Let U be the unitary 
matrix such that C = UQ. Then 

\\iQN-Q)S\\2=\\iUQN-m\2. 

Write UQn = Qn + where Qw is upper triangular and L]y is strictly lower 
triangular, then UQn = Qn + In where qn = HivQwLljY and In = L]^ + Qn — Qn- 
Theorem 3.2 shows that H/w^lb tends to zero exponentially fast. Furthermore since 
Qiv is symmetric, 

Pn =Ql = Q^nQn = {UQn)\UQn) 

={qn + In)\qn + In) 

—Q^n^n P T Q^n^n + I^n^n- 


That is, we have 

Cj^CN Qn^^ ~ ^N^^ Qn^^ In^N- 
Since Qat is uniformly bounded so is qn and we find, 

\\{CI,Cn - qUN)S\\2 < . 

Restricting the above quantities to their finite matrices we note because of the strict 
positivity of P, UCArlh is uniformly bounded from below hence is uniformly 
bounded. Furthermore since Cn has the same size as qn, 

ll(^ - < A'sA'". 

where Sn = CnS for any d with finitely many nonzero entries. Note that the 
factor QnC^^ is upper triangular while is lower triangular. The above 

inequality shows that {{{qnC^^ — /)h ||2 < KsX^. This completes the proof. □ 


§4. Numerical Examples 

In this section we give three examples to illustrate how the computational 
method works for polynomial matrix factorizations. 

Example 4.1. We first consider a univariate polynomial matrix 

8 + 2 + 1/^ 1 + z' 

1 + lA 1 


P{z) := 


12 



It is clear that the matrix is Hermitian and positive dehnite. We write 


Piz) 


8 

1 


1 

1 



2 ; + 


1 

1 


0 

0 


/z. 


We assemble a bi-inhnite Toeplitz matrix whose 10 x 10 block is as shown below. 


-81110 
110 0 0 
10 8 11 

10 110 

0 0 10 8 

0 0 10 1 

0 0 0 0 1 

0 0 0 0 1 

0 0 0 0 0 

.0 0 0 0 0 


0 0 0 0 0 - 

0 0 0 0 0 

1 0 0 0 0 

0 0 0 0 0 

1110 0 
1 0 0 0 0 

0 8 111 

0 110 0 

0 10 8 1 
0 10 11. 


We nse the Cholesky decomposition method to a 20 x 20 central block and get a 
lower triangnlar matrix F. Let PO be the right and bottom 2x2 block from F 
which is 

■ \/385 Q 

pn .= 7 _ 

^ ^ • 6 \/2310 ■ 

. V385 55 , 

Choose the 2x2 block next to PO as follows 


PI : = 


\/385 -\/2m0 

55 3 85 

\/385 -\/2m0 

. 55 385 


Dehne Q^{z) = PO + Pl/z and then we have P{z) = Q{zyQ{z). □ 

Example 4.2. We next consider a bivariate polynomial 

P{x, y) = 41 + 5x^ + 5y^ + 15/a; + 20/y + 5/x^ + 5/y^ + 15a: + 20y + 5xy 

+ 8y/x + 5/{xy) + 8x/y + 2xjy^ + Oyjx^ + 8x^ jy + x^ jy^ + 2y^ jx + y^ jx^ 


It is a positive polynomial since P{x,y) — p{x,y)p{l/x,l/y) with p{x,y) = 5 + 
2x + 8y + xy + x‘^ + y'^. Let ns write 


Pix,y) = [l,l/x,l/xyP{y) 


X 

2 

X 


with 


P(.y) ■= 


il 

3 


+ 


5y 

3 

15 


, 20 , _5_ 
3y^ 3y2 
5 


+ Wy 


^ + 4y+^+y^ 
5 + 3y + y^ 


f + ly + t + ¥ 


3 

15 

2 


3y 
+ 4y + 


2y 


3y2 


5 + —I —2 

y y-^ 

15 I 5 , 4 , ^ 

2 2^ ^ ^ 7/2 


2 

41 I 5,,2 , 20 


t I O 

3y ^ 37/2 


+ Wy 
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Then we can write 


P{x,y) = [l,l/x,l/x^]P{y) 


X 


where P{y) = Yl‘j=- 2 Pjy^ with p_ 2 , ■ ■ ■ ,P 2 being given below: 


Po = 


P2 = 


« 15 r 

k ik ^ 

^ i A 

^23 


■ 

■ 20 

5 

0 ' 

,P 1 = 

3 

4 
_ 3 

1, 

3 

4 

5 

1) 

3 . 


,P-1 =Pl, 


0 0 


,P-2 =pI- 


We now assemble a bi-inhnite Toeplitz matrix whose 9x9 central block are shown 
as follows: 


- 1 

15 

2 

5 

20 

3 

5 

2 

0 

5 

3 

0 

0 

15 

41 

15 

4 

20 

5 

1 

5 

0 

2 

3 

2 

3 

2 

3 

5 

15 

2 

41 

3 

3 

4 

20 

3 

1 

1 

5 

3 

20 

3 

4 

3 

41 

3 

15 

2 

5 

20 

3 

5 

2 

0 

5 

20 

4 

15 

41 

15 

4 

20 

5 

2 

3 

2 

3 

2 

3 

2 

0 

5 

2 

20 

3 

5 

15 

2 

41 

3 

3 

4 

20 

3 

5 

3 

1 

1 

20 

3 

4 

3 

41 

3 

15 

2 

5 

0 

5 

1 

5 

20 

4 

15 

41 

15 

3 

2 

3 

‘-± 

2 

3 

2 

0 

0 

5 

Q 

0 

5 

0 

20 

Q 

5 

15 

0 

41 

Q 


We nse the Cholesky factorization of a central block matrix of size 60 x 60. Let F 
be the lower triangnlar factorization. Then choose Qo to be the 3x3 block at the 
bottom and right of F, Qi the 3x3 block next to Qi and Q 2 the 3x3 block next 
to Qi that is 


Qo = 


Qi = 


Q2 — 


3.185602126 

1.873651218 

1.524622962 


0 

2.539725049 

1.128505745 


0 

0 

2.269126602 


1.797364251 0.08381502303 -0.0003518239229 

0.7675275947 1.633796832 0.06150315980 

0.00008111923034 0.9665117592 1.856367398 

0.5231873284 0.007768330871 0.08530594055 
0 0.6562390159 0.1143305535 

0 0 0.7344969935 


Let ( 5 (y)t = QQ + Qily + Q 2 ly'^ and then Qiy^Qiy) ^ P{y)- In fact the maximnm 
error of each entry of — P{y) is less than or eqnal to 10 “®. □ 
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Example 4.3. Let us consider a bivariate polynomial which has a zero on the 
bi-torus: 


P{x, y) = 30 •+- 14/a: + 11/y + ^jxjy -|- 14a: -|- ^xjy -|- lly -|- 6y/x -|- Axy. 

It is the product of P{x, y) — (4 -|- 3a: -|- 2a/ -|- 1)(4 -|- 3/a: -|- 2/a/ -|- 1) which is zero at 
X = —1, y — —1. We write 


P{x,y) = po{y) +pi{y)x + p-i{y)/x 


for po(y) = 30 -f- ll/a/ 11a/, pi{y) = IA + Qy + 4/a/, and p-i{y) = IA + Ay + 6/a/. 
It is the symbol of an bi-inhnite Toeplitz matrix. One of its central section is as 
shown below. 


-ll/a/+ 30 + 11a/ 
6 /a/ + 14 + 4a/ 

0 

0 


4/a/ + 14 + 6a/ 
ll/a/+ 30 + 11a/ 
6 /a/ + 14 + 4a/ 

0 


0 

4/a/ + 14 + 6a/ 
ll/a/+ 30 + 11a/ 
6 /a/ + 14 + 4a/ 


0 

0 

4/a/ + 14 + 6a/ 
ll/a/+ 30 + 11a/. 


Since P{x, y) has no simple factors (see the next section), any central sections of the 
bi-inhnite Toeplitz matrix is positive by Lemma 5.1. We consider several central 
sections Pm of size an = 16 x 16, 32 x 32, 64 x 64 and 128 x 128. For each of 
these central sections. Pm is a univariate polynomial in y with matrix coefRcients 
and Pm{y) is positive. Thus, Pm{y) = Qm{yyQm{y)- To compute Qm{y), we use 
the computational method in §3 to yield an approximation Qm of Qm- As the size 
of central sections increases, the Qm converges to the corresponding entries in the 
bi-inhnite Toeplitz matrix. We use the entries on the last row of Qm to construct 
an approximation of Qm{y) and hence the factorization of P{x, y) and listed below. 


size 

16 X 16 
32 X 32 
64 X 64 
.128 X 128 


factorization 

4.01207952 + 2.984741799a: + 2.000226870// + 0.996712925a:// 
4.004041536 + 2.994924757a: + 2.000034879// + 0.998949058a:// 
4.001381387 + 2.998269650a: + 2.000005690// + 0.999648058a:// 
4.00069369 + 2.999134582a: + 1.99999896// + 0.999821915a:// . 


As we know that the factorization is 4 + 3a: + 2// + 1, the approximations are very 
good. □ 


§5. Nonegative bivariate Trignometric Polynomials 

Finally we consider the problem of factorization of nonnegative multivariate poly¬ 
nomials. Let us start with P{z,w) > 0. If for some zq with \zq\ = 1, P{zo,w) = 0 
for all w with |r(;| = 1, we say that P{z, w) has a simple factor at zq. If P{z, w) has 
a simple factor at zq, then P{z,w) has factors {z — zq) and Qjz — 1/zq). Let us 
factor them out. Then P{z^w)/{{z — zq){1/ z — 1/ zq)) is still nonnegative. Similarly, 
if P{z,wo) = 0 for all z with \z\ = 1, P{z,w) has a simple factor at wq. In this 
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case, P{zj w) has two factors {w — wq) and {1/w — I/wq) which can be factored ont 
from P{z, w). Withont loss of generality, we may assnme that P{z, w) >0 does not 

n 

have any simple factors. Writing P{z,w) = pj{w)z^, we view that P{z,w) is 

j = -n 

a polynomial of ^ and it is the symobl of a bi-in£nite Toeplitz matrix in (4) with 
w in place of Z 2 - We have the following 

Lemma 5.1. Suppose that P{z, w) > 0 does not have any simple factors. Then any 
central section of the bi-inhnite Toeplitz matrix in (4) is strictly positive definite. 

Proof: Since P{z, w) > 0, we know that any central section of the matrix in (4) 
is nonnegative dehnite. Snppose that a central section Tm{w) of the matrix in 
(4) is not positive definite for w = wq. Then there exists a vector x snch that 
Tra{wo)yi = 0, i.e., x^T^(m;o)x = 0. Thns, we have, for z 

1 

0 = x^T^(m;o)x = — F{x)*P{z, wo)F{x)d9. 

ZTT Jq 


It follows that 

\F{x)\‘^P{z,wo) = 0, a.e. 

and hence, P{z^ wq) = 0 since |T'(x) | ^ 0, a.e. and P(z, wq) is a Lanrent polynomial. 
That is, P{z,w) has a simple factor at wq. This contradicts the assnmption on 
P{z,w). □ 

Thns, for a central section Pm of size m x m in the matrix in (4), Pm is positive. 
Since Pm is a matrix polynomial in rc, by the matrix Fejer-Riesz factorization theo¬ 
rem (cf. [Helson’64]), Pm can be factorized into Qm, he., Pmiw) = Qmiw)"^Qmiw). 
Intnitively, the polynomial Qm is a good approximation of the factorization of the 
bi-inhnite Toeplitz matrix V in (4) as m snfhciently large. In the previons section, 
we presented an example (Example 4.3.) of P{z, w) which is nonnegative withont 
simple factors. Using onr symbol approximation method, we compnte an approx¬ 
imation of the factorization of Pm for m = 16, 32, 64, and 128. The nnmerical 
compntation shows the factorizations converge. 

Let ns now discnss the convergence a little bit more in detail. For simplicity, let 
.4. be a bi-inhnite Toeplitz matrix associated with a nnivariate Lanrent polynomial 
A(z) and An = icijk)-N<j,k<N be a central section of size {2N -|- 1) x {2N -|- 1) 
for a positive integer N. Snppose that each is strictly positive. Thns we can 
obtain a factorization = Pn^n by Cholesky factorization. 

Lemma 5.2. For any x,y G £ 2 , xn^^tvY := xj^^Tvyw converges to x^^y as 
N —-l-cxD, where xn — {x-n, ■ ■ ■ , xq, ■ ■ ■ , xat)^ is the central section of size 2N + 1 
of X around the index 0 and similar for hn. 
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Proof: For an integer > 0, 


x^^Ary - 

/•27r 

r.27r 

(F(x7v) - F{yi))* A{z)F{yN)de 


(F(xAr)M(z)F(yAr) - F(x)M(^)F(y)) dO 


'0 


1 


'0 


1 


r*27r 


+ — / F(x)M(z)(F(yA.)-F(y))rf^ 


where z = e*^. In the hrst ineqnality we nsed the fact that x^^^AfX = (nArx)^^^nArx 
where Hat is the projection dehned in section 3. Thns 


|xMwy - xMy| 

<||x - XAr|| 2 P( 2 :)||oo||y ||2 + ||y - yw|| 2 P(^)||oo||x ||2 
-^0 


as N ^ + 00 . Here, ||H(2)||oo denotes the maximnm norm of A(z) over the circle 
1^1 = 1. This completes the proof. □ 

A conseqnence of the above Lemma 5.2 is that ||i3Arx||| converges to x^Ax. If 
A can be factored to A = Then ||i3A/^x||2 —^ ||i3x||2. The following is another 
conseqnence of Lemma 5.2. 

Lemma 5.3. Let Bn be a factorization of An, be., An = B]^Bn- Then \\Bn\\ is 
bounded independent of N. 

Proof: By Lemma 5.2, there exists a constant C snch that for N large enongh, 

||HArx ||2 = x^Aatx < x^Ax + C = ||x|| 2 ||A( 2 :)||oo + C. 

Hence, HHatH := max ||i3A/'x||2 is bonnded. □ 

Il^ll2=l 

Note that all Bn are banded npper triangnlar matrices with one half the band 
width as that of A. Thns, each row (or colnmn) of Bn has hnitely many nonzero 
entries. Lemma 5.3 implies that each row (or colnmn) of Hat is bonnded in £2 norm 
and hence each entry in any row is bonnded. Therefore there exists a snbseqnence 
of Bnj snch that each entry with indices (j, k) in BNi converges as i —+cx). That 
is, for any vector x = {xi)i ^2 £ £2 with hnitely many nonzero entries XiS, we have 

Hat^x —> Hx. 

for a bi-inhnite matrix B. By Lemma 5.2 again, we have x^H^Hy = x^Ay for all 
vectors x and y with hnitely many nonzero entries. However since these are dense 
in £2 we hnd B^B = A. Note that B is an npper triangnlar matrix with the half 
the band width as that of A. If H is a Toeplitz matrix, we immediately know that 
A{z) has a factorization snch that A{z) = B{z)*B{z). Therefore, we end with 
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Theorem 5.4. Let P{z,w) be a nonnegative Laurent polynomial with no simple 
zeros. Let V be a bi-infinite Toeplitz matrix with Laurent polynomial entries in 
w, Pn be the central section as described above and be its upper triangular 
Cholesky factor. Then there is a subsequence of Bn convergent to B entrywise, 
where V = B'^B. If B is Block Toeplitz, then P{z, w) can be factored into a sum of 
square magnitudes of finitely many polynomials in z and w. 

Theorem 5.4 provides a computational method to check if a nonnegative Lau¬ 
rent polynomial P{z, w) can be factorized. That is, we compute Cholesky factoriza¬ 
tion of central sections of the bi-infinite Toeplitz matrix V associated with P{z, w) 
and observe if the factorization matrices converge to a Toeplitz matrix or not. If 
they converge, P{z,w) has a factorization. 


§6. Remarks 

1. It is interesting to point out that the symbol approximation method discussed 
in [Lai’94] is very much like the Bauer method invented in 1955 (see [Sayed 
and Kailath’Ol] and its references). A slight difference is that the singular 
value decomposition (SVD) instead of the Cholesky decomposition is used to 
factorize the matrices associated with Laurent polynomial P{z) > 0. 

2. When P{z) is a matrix polynomial in the univariate setting [Hardin, Hogan and 
Sun’04] have demonstrated a constructive method to factor P{z) = Q^zYQ^z) 
when P{z) has a nonzero monomial determinant. 

3. When P{z) is a matrix polynomial in the univariate setting [Youla and Kazan- 
jian’78] used a Bauer type method to compute the factorization of P{z). They 
showed that the solution obtained from the Bauer type method converges to 
the factorization under a weaker condition that 


1 


logdet P{z)d9 > —oo 


than the positivity condition of P{z). In our Theorem 3.1. we showed the ex¬ 
ponential convergence of the method which greatly improves their convergence 
analysis. 
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